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Abstract Baryons made of two or three heavy quarks can be described in the modern 
language of non-relativistic effective field theories. These, besides allowing a rigorous 
treatment of the systems, provide new insight in the nature of the three-body interac- 
, tion in QCD. 
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1 Motivations 
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Baryons made of two or three heavy quarks offer an interesting alternative to quarko- 
nium for studying the dynamics of non-relativistic systems in QCD and for investigating 
the transition region from Coulombic to confined bound states. 

The modern approach to quarkonium physics consists in taking advantage of the 
hierarchy of non-relativistic energy scales in the system by constructing a suitable hi- 
erarchy of effective field theories (EFTs) pQ. The energy scales are the heavy-quark 
mass, m, the typical momentum transfer, p <C m, whose inverse sets the typical dis- 
tance, r, between the heavy quark and the antiquark, and the typical kinetic energy, 
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, E <C p, whose inverse sets the typical time scale of the bound state. In the ultimate 

■^j- ■ EFT, obtained after integrating out gluons of energy and momentum of the order of 

m and p, the interaction between heavy quarks is organized as an expansion in powers 
of I/m and r. At zeroth-order in r, the interaction is entirely encoded in the quark- 
antiquark potential, which, at zeroth-order in I/m, reduces to the static potential. 
Terms proportional to powers of 1/m and r can be systematically added. 

An analogous hierarchy of energy scales characterizes also baryons that contain 
at least two heavy quarks. Hence, we may describe these systems by means of EFTs 
analogous to the ones suited for heavy quarkonium [2][3]. In particular, the ultimate 
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EFT, obtained after integrating out gluons of energy and momentum of the order of 
the heavy-quark masses and of the typical momentum transfer between heavy quarks, 
is organized as an expansion in the inverse of the heavy-quark masses and in the 
distance between the heavy quarks. For equal heavy-quark masses, the structure of the 
Lagrangian and its power counting are similar to the quarkonium case; the Lagrangian 
reads 

C = -\F^F a " u + J2 If iJpqf + £ 5C^ , (1) 

/ i,j 

where 8C^ 1 '^ are terms containing the heavy quark or antiquark fields, which are 
proportional to 1/m 1 x [typical distance between heavy quarks/antiquarks]-' . The fields 
qf are nj light-quark fields, assumed to be massless. The heavy quark or antiquark 
fields in 5£y^' are twice those in the bound state, e.g, in the quarkonium case, 5£^'^ 
contains two quark and two antiquark fields. For different heavy-quark masses, more 
scales are involved; in the following, we will not consider such cases. 

Although the EFTs for quarkonium and baryons made of two or three heavy quarks 
are similar in structure, they are characterized by different degrees of freedom. This is 
best seen when the typical distance between heavy quarks is smaller than the inverse 
of the hadronic scale Aqqd, which is the case that we will discuss in this note. At these 
distances, gluons may resolve coloured degrees of freedom. These are gluons or light 
quarks or, in the quarkonium case, QQ states in a colour singlet or in a colour octet 
configuration. In the case of baryons made of two heavy quarks Q and a light quark 
q, at distances smaller than 1/^4qcd> gluons can resolve gluons, light quarks and QQ 
pairs in a colour antitriplet or in a colour sextet configuration. It is the binding of the 
antitriplet with the light quark q that forms the QQq baryon. The system very much 
resembles a heavy-light meson, with the heavy antiquark replaced by a QQ antitriplet; 
this fact may be exploited to deduce some properties of the QQq baryons from the 
corresponding Qq mesons [4]. In the case of baryons made of three heavy quarks Q, 
at distances smaller than 1/Aqq^,, gluons can resolve gluons, light quarks and QQQ 
states either in a colour singlet or in two different colour octets or in a colour decuplet 
configuration. 

The lattice evaluation of the QQQ static potential has a long tradition (see e.g. [5] 
and references therein), while the static potential between a QQ pair in the presence of 
a light quark has been evaluated on the lattice only recently [6] . Expressions for the 1/m 
and the 1/m 2 spin-dependent QQQ potentials in terms of Wilson loops can be found 
in [2J, but have not been calculated on the lattice yet (while the complete expressions of 
all the QQ potentials up to order 1/m 2 can be found in [7] and the most recent lattice 
determinations are in [8]). Perturbative studies of the static potential may help to un- 
derstand the transition region from the perturbative to the non-perturbative regime. In 
the case of the QQ potential, this region is characterized by the smooth transition from 
a Coulomb potential to a linear raising one (for recent studies, see [S]). In the case of 
the QQQ static potential, the transition from the perturbative to the non-perturbative 
regime is accompanied by the emergence of a three-body potential that depends on 
one length only (see, for instance, [5]). This is a rather spectacular phenomenon, which 
has been investigated only recently from a perturbative perspective [10J . 

There is so far no experimental evidence of QQQ baryons, while a few years ago the 
SELEX experiment has claimed evidence of possible doubly charmed baryon states [11] . 
Until today this evidence has not been confirmed by other experiments, but is mostly 
behind the revival of interest in this kind of systems at the mid of this decade. 
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2 EFT for QQq 

The EFT Lagrangian that describes QQq baryons below the momentum transfer scale 
p, assumed to be larger than /IqcDi has the general form of Eq. {TJ, with SC^'^ made 
out of four quark fields [2]. 
The term <5£ (0 ' 0) is 

6C^ = J d 3 r Tt [iD - V± 0) ] T + [iD - V^] E, (2) 

where T = (T 1 ,T 2 ,T 3 ) are the three independent QQ antitriplet fields, E = (E 1 , E 2 , 
. . . , E e ) are the six independent QQ sextet fields and the gauge fields in the covariant 
derivatives acting on the antitriplet and sextet fields are understood in the 3 and 6 
representations respectively. The matching coefficients and can be identi- 

fied with the antitriplet and sextet static potentials respectively. The QQ antitriplet 
static potential has been calculated recently up to next-to-next-to-leading order 
(NNLO) [ID]; it reads 

v ( )_ 2 a s (l/H) f 1 , q s (1/H) [31 , 00 _ fW , 4_ 
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+ Iy 7b + 27 7b + 27" + ^rH)' (3) 

where a s is the strong-coupling constant in the MS scheme. The QQ sextet static 
potential is repulsive at leading order: V^P = a s /(3|r|). 

Several terms contribute to . The one responsible for the hyperfine splitting is 

5£ (1,0) = ]4^BT r t£F (V 1) + a (2)\ . gB a T a ^ (4) 

where vfi£) BT = 1 + 0(a 2 ) is a matching coefficient of the EFT, cp = 1 + . . . is the 
matching coefficient of the chromomagnetic interaction in the heavy quark effective 
theory, which is known up to three loops [12], cr^' is a Pauli matrix acting on the 
heavy quark labeled i and are the Gell-Mann matrices in the 3 representation. 
From Eq. Q, it follows that the hyperfine splitting between the S-wave ground state 
of a doubly heavy baryon of spin 1/2 {Sqq) and the corresponding state of spin 3/2 
{Sqq) may be related to the hyperfine splitting between the S-wave ground state of a 
heavy-light meson of spin (Pq>) and the corresponding state of spin 1 (P<g/): 



- Mb 



Q (Q) 



QQ ~QQ 4m (Q 



~ { Mp Q> ~ Mp Q' ) 



l + O [a 



2 ^QCD ^QCD 



<Q 



(5) 



where we have kept different the mass, ttiq, of the heavy quarks in the baryons from 
the one, tuqi, in the mesons. The obtained figures compare well with existing lattice 
determinations (see discussion and references in [2], a more recent unquenched lattice 
determination of Mo* — Ms bb , which is consistent with previous quenched determi- 
nations, may be found in |13j). 

Other terms in the effective Lagrangian have been derived in [2]. 
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3 EFT for QQQ 

The EFT Lagrangian that describes QQQ baryons below the momentum transfer scale 
p, assumed to be larger than vIqcDi nas the general form of Eq. ([1]), with SC^'^ made 
out of six quark fields [2] . 
The term <5£ (0 ' 0) is 
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where S is the singlet field, O = {^) s ^j ' With ° A = (0 A1 , A2 , . . . ,0 Aa ) and 

O = (O , O , . . . , O ), are the fields that parameterize the two possible octet 
configurations of three quarks, A = (A 1 , A 2 , . . . , A 10 ) are the ten independent QQQ 
decuplet fields and the gauge fields in the covariant derivatives acting on the octets 
and decuplet fields are understood in the 8 and 10 representations respectively. The 
matching coefficients Vg°\ Vq^ and vj^ can be identified with the singlet, octet and 
decuplet static potentials respectively; note that the octet potential Vq^ is a 2 x 2 ma- 
trix. The QQQ singlet static potential has been calculated up to NNLO in |1U| ; 
it reads 
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+y| body (r 1 ,r 2 ,r 3 ), (7) 

where n, r2 and r 3 are the distances between the heavy quarks; only two of them are 
independent: if we call x\, x 2 and a; 3 the coordinates of the three quarks, our choice 
is r\ — xi — X2, r 2 = x\ — a; 3 and r 3 = x 2 — a; 3 , which implies that r 3 = r% — r\. 
yS o y . g ^ e three-body part of the perturbative potential, defined as the part of the 
potential that vanishes when putting one of the quarks at infinite distance from the 
other two. 

The three-body part is a specific feature of the two-loop potential; it reads 
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Fig. 1 The normalized three-body potential, 2Vg body /ctg, plotted as function of L in arbi- 
trary units for the two geometries described in the text. 



with R — xp — y\, R — \R\ and M = \p\ \J x(l — x) + \X\^Jy(l — y). This term has a 
different dependence on the positions of the three quarks with respect to the Coulomb 
potential. It is finite for all configurations: it vanishes when one of the quarks is pulled 
at infinite distance (note that in this case vjP^ becomes V^ 0) ) and still remains finite 

when two quarks are put in the same position (note that in this case Vg ' becomes the 
QQ singlet static potential). Hence, its dependence on the geometry is much smoother 
than the dependence on the geometry of the Coulomb potential, although it still clearly 
depends on it. These features can be seen by plotting vJ body as a function of L for 
two different set of configurations; L is the sum of the distances of the three quarks 
from the so-called Fermat (or Torricelli) point, which has minimum distance from the 
quarks. Figure [JJ a) shows the three-body potential when we place the three quarks in 
a plane (x, y), fix the position of the first quark in (0, 0), the position of the second one 
in (1,0) and move the third one in (0.5 + 0.125 n x , 0.125 n y ) for n x € {0, 1, . . . , 20} and 
n y £ {0, 1, . . . , 24}. The plot shows a clear dependence on the geometry at fixed L (i.e. 
different configurations with same L give different potentials), although weaker than 
in the Coulombic case. Figure [Tib) shows the three-body potential in a geometry used 
in [5] , which consists in placing the three quarks along the axes of a three-dimensional 
lattice: (n x , 0, 0), (0, n y , 0) and (0, 0, n z ), for % £ {0, 1, ... , 6} and n y ,n z € {1, ... , 6}. 
The plot shows a weaker dependence on the geometry than in the previous case. As 
the comparison with Fig. (Ha) indicates, the weaker dependence is an artifact of the 
chosen configurations rather than a physical effect. Clearly, the precise identification of 
the transition region from a two-body dominated potential to a three-body dominated 
one would require to perform lattice calculations in geometries different from the one 
used in Fig.[ljb). The use of different geometries could also be important to assess the 
nature of the long-range three-body potential. 

The one-gluon exchange mixes the octet fields, so that the octet potential Vq^ 
is a non-diagonal 2x2 matrix already at leading order. Choosing O and O to be 
respectively symmetric and antisymmetric for exchanges of the quarks located in ii 
and X2 (a different choice of the octet fields would correspond to a field redefinition 
leading to a different octet potential), we obtain 
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while at leading order the decuplet potential is 



A 



(0) _ «s 




) 



(11) 



We mention that there exist lattice data that show very clearly the singlet, octet and 
decuplet QQQ static potentials, although in an equilateral geometry, where the two 
octets are degenerate [14] . 

Other terms in the effective Lagrangian have been derived in [2J. 

4 Conclusions 

Systems made of two or three heavy quarks or antiquarks develop similar hierarchies 
of energy scales and may be treated in similar EFT frameworks. 

In the case of QQ mesons, the static potential has been determined up to next- 
to-next-to-next-to-leading order in perturbation theory and to a high accuracy on the 
lattice. The Coulomb behaviour starts getting substantial modifications at distances 
around 0.2 fm turning over a linearly raising potential at larger distances. Therefore 
the transition region can be studied to a large extent with perturbative methods. 
Terms proportional to powers of 1/m and r have been also calculated and included 
systematically in physical observables. 

In the case of QQq baryons, the static potential has been determined up to NNLO 
in perturbation theory and recently also on the lattice. Terms proportional to powers 
of 1/m and r in the Lagrangian have been matched (mostly) at leading order and used 
to determine, for instance, the expected hyperfine splitting of the ground state of these 
systems. If early experimental evidences will get confirmation in future experimental 
facilities, this will mark the beginning of a future new spectroscopy. 

Finally, in the case of QQQ baryons, the static potential has been determined 
up to NNLO in perturbation theory and also on the lattice. The transition region 
from a Coulomb to a linearly raising potential is characterized in this case also by 
the emergence of a three-body potential apparently parameterized by only one length. 
While we have argued that more lattice studies employing different geometries would 
be necessary to precisely identify the transition region, we have also shown that in 
perturbation theory a smooth genuine three-body potential shows up at two loops. 
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